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We present parameter-free predictions of the nuclear modification factor, Raa(Pt, s), of high pr 
pious produced in Pb+Pb coUisions at ^/s^J^ = 2.76 and 5.5 ATeV based on the WHDG/DGLV 
(radiative+elastic+geometric fluctuation) jet energy loss model. The initial quark gluon plasma 
(QGP) density at LHC is constrained from a rigorous statistical analysis of PHENIX/RHIC vr" 
quenching data at \/sjv]v ~ 0-2 ATeV and the charged particle multiplicity at ALICE/LHC at 
2.76 ATeV. Our perturbative QCD tomographic theory predicts significant differences between jet 
quenching at RHIC and LHC energies, which are qualitatively consistent with the pT-dependence 
and normalization — within the large systematic uncertainty — of the first charged hadron nuclear 
modification factor, R'aaj data measured by ALICE. However, our constrained prediction of 
the central to peripheral pion modification, Rcpipr), for which large systematic uncertainties 
associated with unmeasured p+p reference data cancel, is found to be over-quenched relative to 
the charged hadron ALICE _R^p data in the range 5 < pt < 20 GeV/c. The discrepancy challenges 
the two most basic jet tomographic assumptions: (1) that the energy loss scales linearly with the 
initial local comoving QGP density, po, and (2) that po oc dN'^^ {s,C)/dy is proportional to the 
observed global charged particle multiplicity per unit rapidity as a function of ^/s and centrality 
class, C. Future LHC identified {h — Tv,K,p) hadron Raa data (together with precise p-l-p, 
p-|-Pb, and Z boson and direct photon Pb-|-Pb control data) are needed to assess if the QGP pro- 
duced at LHC is indeed less opaque to jets than predicted by constrained extrapolations from RHIC. 
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I. INTRODUCTION 

The first LHC Pb+Pb data at 2.76 ATcV [1- 

6] provide important new consistency tests of dynami- 
cal models developed over the past two decades to pre- 
dict multiparticle observables in nuclear reactions at 
^/s = 0.2 AGeV from the Relativistic Heavy Ion Col- 
lider (RHIC/BNL) [7-10]. The new LHC energy fron- 
tier with Pb+Pb probes the physics of strongly coupled 
Quark Gluon Plasma (sQGP) [11-13] at densities approx- 
imately twice as high as at RHIC, with temperatures up 
to T ~ 500 MeV, well above the deconfinement transition 
region predicted by lattice QCD [14, 15]. 

The theoretical understanding of the dynamical prop- 
erties of this new form of matter is however far from 
complete. On the one hand, bulk radial and differential 
azimuthal anisotropic elliptic flow data of low transverse 
momenta {p ~ STq sa 1 GeV) partons were found to 
be consistent with near "perfect fluid" flow and suggest 
highly nonpcrturbative physics of this new form of matter 



[16, 17]. This has led to proposals that the bulk proper- 
ties of the sQGP may be better approximated via strong 
coupling (supergravity dual) holography models [18-29] 
than via perturbative QCD (pQCD) based quark and 
gluon quasiparticle Hard Thermal Loop (HTL) [30, 31] 
approximations. 

On the other hand, short wavelength {pT ^ 10 — 20 
GeV/c) properties of the sQGP, as measured via the nu- 
clear modification factor of pions [32-34], were found to 
be well predicted by pQCD based HTL partonic radiative 
energy loss theory [35-57]. The nuclear modification fac- 
tor of high transverse momentum hadron, h, fragments 
ill A + B h + X and centrality class C used to probe 
the short wavelength dynamics in an sQGP is defined as 



RAB{y,PT;Vs,C) = 



dN^+^^f'{y,PT, ^s,C)ldy<fpT 
TAB{C)daP+P^'^{y^)/dyd^PT 
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(1) 

For a fixed \/s center of mass (cm) energy (per nucleon 
pair) and nucleon-nucleon (NN) inelastic cross section 
'^n'n(Vs) the mean number of elementary binary NN col- 
lisions in centrality class C is given by f7™j^TAA, where 
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[58] 



Tab{C) = /J(fx±TA{xA 



^)TBix^ + l)) (2) 
/ bee 



in terms of the Glauber nuclear thickness profile 
TAix±) ~ J dzpA{z,xj_) and Wood-Saxon nuclear den- 
sity Pa normalized to A. 
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FIG. 1. WHDG model [53] predictions (blue bands extrapo- 
lated from the RHIC constrained green band) for the nuclear 
modification factor of 7r° in Pb-f Pb 2.76 ATeV LHC are com- 
pared to ALICE/LHC [1] c/tanjed hadron nuclear modification 
data in central (red solid) and peripheral (open red) reactions. 
The PHENIX/RHIC Au+Au-> 7r° nuclear modification data 
[34] are shown by black dots. The brown triangles and blue 
stars represent the charged hadron PHENIX [32] and STAR 
[33] data, respectively. The blue band of WHDG predictions 
corresponds to the l-cr medium constraint set by PHENIX [34] 
extrapolated to LHC via the ALICE charged particle rapidity 
density [2]. The wide yellow band is the current systematic 
error band of the (red dot) LHC data due to the unmeasured 
p+p reference denominator. 

In the absence of both initial state and final state nu- 
clear interactions Rab = 1- For pT below some charac- 
teristic medium dependent transverse momentum "sat- 
uration" scale, Qs{pt, \/si A), the initial nuclear par- 
tonic distributions functions (PDFs) [59-61] fa/A{x — 
'^PT/y/s,Q^ ^ Pt) < ^fa/N{x,Q^) are expected to be 
shadowed, leading to Raa < 1 because the incident flux 
of partons is less than A times the free nucleon parton 
flux. Color Glass Condensate (CGC) models [11, 62- 
68] have been developed to predict Qs{pt, -\/s, A) related 
initial state effects from first principles. While the mag- 
nitude of Qs at LHC is uncertain and will require future 
dedicated p-|-Pb control measurements to map out, cur- 
rent expectations are that < 5 GeV at LHC in the 
central rapidity region. This should leave a wide jet to- 



mographic kinematic window 10 < px < 200 GeV in 
which nuclear modification should be dominated by final 
state parton energy loss and broadening effects. In this 
paper, we therefore assume that initial state nuclear ef- 
fects can be neglected in the 10 < px < 20 (i.e. x > 0.01) 
range explored by the first ALICE data [1] . We note that 
from Fig. 1, and as discussed in detail below, our RHIC 
constrained jet quenching due to final state interactions 
alone already tends to over-predict the pion quenching 
at LHC and therefore leaves no room for large addi- 
tional shadowing/saturation effects in the [68-70] in this 
> 100 GeV^ kinematic window — unless the sQGP is 
much more transparent at LHC than expected from most 
extrapolations of jet quenching phenomena from SPS and 
RHIC to LHC energies. 

The main challenge to pQCD multiple collision theory 
of jet tomography and AdS/CFT jet holography is how to 
construct a consistent approximate framework that can 
account simultaneously for the beam energy dependence 
from SPS to LHC energy and for the nuclear system size, 
momentum, and ccntrality dependence from p-\-p to U + 
U of four major classes of hard probe observables: (1) the 
light quark and gluon leading jet quenching pattern as a 
function of the resolution scale pt, (2) the heavy quark 
flavor dependence of jet flavor tagged observables, and (3) 
the azimuthal dependence of high px particles relative to 
the bulk reaction plane determined from \ow-pT elliptic 
flow and higher azimuthal flow moments, Vnipr); and (4) 
corresponding di-jet observables. 

The first LHC heavy ion data on high transverse mo- 
mentum spectra provide an important milestone because 
they test for the first time the density or opacity depen- 
dence of light quark and gluon jet quenching theory in a 
parton density range approximately twice as large as that 
studied at RHIC. The surprise from LHC is the relatively 
small difference observed between the RHIC [32-34] and 
ALICE [1] LHC data on Raa{10 < pr < 20 GeV), as 
shown in Fig. 1. In addition, there is little difference 
from RHIC to LHC between the differential elliptic flow 
probe, V2{pt < 2), as reported in [3]. The rather striking 
similarities between bulk and hard observables at RHIC 
and LHC pose significant consistency challenges for both 
initial state production and dynamical modeling of the 
sQGP phase of matter. 

In this paper, we focus on the puzzle posed by the 
similarity of inclusive light quark/gluon jet quenching at 
RHIC and LHC by performing a constrained extrapola- 
tion from RHIC using the WHDG model [53] to predict 
at 2.76 ATeV cm energy. We update our earlier 

2007 LHC predictions in [71, 72], by extrapolating the 

2008 l-cr PHENIX/RHIC constraints [34] of the opac- 
ity range at = 0.2 ATeV using the new 2.76 ATeV 
ALICE/LHC [2, 4] charged hadron rapidity density data, 
dNch/drj = 1601±60, in the 0—5% most central collisions 
and 35 ± 2 in the 70 — 80% peripheral collisions. 

We note that in strong coupling AdS / CFT approaches 
to hard jet probes, the pQCD high-p^ jet tomogra- 
phy theory is replaced by a gravity dual jet holographic 
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model. That approach is based on the assumption that 
the 't Hooft couphng, A = AnasNc, as weU as are 
large enough that an approximate lOD supergravity dual 
description of the dynamics may be used. Jet quench- 
ing is then mapped into the problem of classical string 
drag over a black brane horizon in an AdS curved space- 
time with an isometry group that is identical to the ex- 
act conformal symmetry group of the M ~ 4: Supersym- 
metric Yang-Mills (SYM) cousin of QCD. The thermally 
equilibrated strongly coupled supersymmetric QGP is as- 
sumed dual to the black brane. In [73], it was shown that 
with A - 20 and = 3 « 0.5) AdS/CFT holog- 
raphy provides a remarkably robust analytic account of 
both the single nonphotonic electron (heavy quark jet) 
quenching as well as bulk elliptic flow data, a simulta- 
neous description of which has remained out of reach 
of pcrturbativc QCD methods. However, the theoreti- 
cal consistency of heavy quark jet holography remains 
controversial in the literature (see, e.g., [74]). Light jet 
holography of i?^'^ is even more challenging. So far, only 
schematic falling strings scenarios have been proposed to 
address light quark/gluon quenching observables [26, 27]. 
However, high-p^ elliptic moment V2{pt) phenomenol- 
ogy [75-78] appears to require nonlinear path length de- 
pendences, L"-^, more suggestive of falling string sce- 
narios than L"-^ path dependences predicted by pQCD 
for static plasmas. Future heavy flavor tagged jet ob- 
servables [73, 79-81] will help to discriminate between 
competing jet holography vs. pcrturbativc tomography 
models of jet attenuation. 



II. JET TOMOGRAPHY: QUALITATIVE 

One feature common to both pQCD tomography and 
gravity dual holography is that both predict the energy 
loss per unit length, dE/dL, increases monotonically with 
the plasma density or temperature. In this section we 
consider a generic analytic energy loss that can inter- 
polate between a wide range of dynamical scenarios and 
provide qualitative understanding of the quantitative nu- 
merical WHDG results presented in Section III. 

Consider the following power law model for jet energy 
loss [78] 



dP 

^ = -KPV''T2-'^+^(x(r),r) 
dr 



(3) 



where P{t) is the momentum (energy) of a massless jet 
passing through a plasma with a local temperature held 
T{x,t). The solution for an initial energy P(0) and jet 
path, a;(r), a uniform static "plasma brick" of thickness 
L is 



P{0) 



= 1 - KiLT) 



l+b 



T 

m 



(4) 



Wc note that dP{Q)/dP{L) = {P{Q)/P{L)Y is the Jaco- 
bian of the transformation between P(0) and P{L). The 



parameters a and b control the jet energy and path length 
dependence energy (momentum) loss per unit length 
(time) and fixes the power of the temperature or parton 
density p rx dependence. 

The thermal stopping distance Lt, as deflned by 
P(Lt) = Pt = T, is then: 



LT{Pa,T,a,b,c) = 



(1 - a)Kr2-«+'' 



(5) 



In pQCD, the opacity series WHDG/DGLV [42, 43, 47, 
53] for a massless parton jet leads (in a static uniform 
plasma) to dEoLv/dL « -kT^L^ log{E/T). Therefore 
6=1 (or 0) in the LPM regime in static (or Bjorken 
expanding) plasmas. Because E^^^ / log{E) « 1 ± 0.1 in 
the range 5 < E/T < 200, we can simulate pQCD light 
quark jet energy dependence with (a « 1/3, 6=1). The 
density dependence is then roughly linear with T^/^ ^ 
p8/9 jj-^ ^Yie static case. 

Another interesting limiting case in pQCD is for thick 
plasmas, where the deep LPM regime leads to the BDMS 
[45] formula, dEeDMs/dL = -kE^/^L^T^/^ that grows 
with density as p^/^. Increasing the density by a factor 
equal to the ratio of charged particle multiplicities, « 2.1, 
from RHIC to LHC approximately doubles both the GLV 
and BDMS energy loss. 

In distinction to pQCD, the stopping distance for light 
quark jets in the falling string holographic scenario was 
found in [29] to be bounded by Lt < kE^/^/T-^/^. This 
can be simulated via Eq. (5) by choosing a > (2 — 6)/3 
in the range a S [0, 1/3]. The special case 6 = 2 and a = 
1/3 is the one favored phenomenologically [76-78, 82] by 
the high Pt azimuthal anisotropy V2{pt) data. For that 
case, dEAds/dL = -kE^^^L'^T^^^^ . Thus the energy 
loss grows even faster with density, cx p^^^^. We see that 
quite generally dE/dL should increases significantly with 
density. 

For a more realistic Bjorken longitudinally expanding 
plasma of transverse thickness L the density decreases 
approximately as T^{t) cx {dNch/dy)/{TL'^). In this 
Bjorken brick case 



P{L) _ ( ,{dN/dyY^-''+'>y^ 

— - II AC 



P(0) 



(LP(0)) 



l-a 



(6) 



To estimate the variation of Raa in this case consider 
the distribution of initial jet transverse momenta (at mid 
rapidity) approximated as 



da{pT)/dpT oc l/pj(^-'^) 



(7) 



with n{pT, s) = — dlog (d(T(s)) /d\og{pT) given by the cm 
energy dependent parton spectral index computed from 
pQCD. 

The nuclear modification factor for jet partons (ne- 
glecting hadron fragmentation) is simply a change of vari- 
ables which leads with Eq. (6) for a fixed centrality class 
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to 



RAAiPf] s,A) 



dpo da{po{pf))/dp 
dpf d(T{pf)/dp 

,(dA^/dy)(2-«+'')/3 



1 + k' 



[LpfY- 



.(8) 



Given a model specified by (a, b) the coupUng parameter 
K'(a, h) can be fixed by fitting RHIC Raa at one reference 
momentum, e.g. pr ~ 10 GcV/c. 

We show in Fig. 2 the evolution of this simple an- 
alytic model of Raa from RHIC to LHC for the case 
a = 1/3 energy loss and for different 6 = 0, 1, 2 that cor- 
respond to path length scaling for elastic, inelastic and 
AdS falling strings scenarios. For this qualitative plot, 
we take spectral indices (see Fig. 3 of the next section) 
as nnHic = 't-i(pt) ~ 5.5 + {px — 10)/10 that rises by 
1 unit in the pr = 10 - 20 range at RHIC 0.2 ATeV, 
while riLHC = f^2 ~ 4.5 approximately independent of 
PT at LHC 2.76 ATeV. We fix n' for each b by demand- 
ing RAAipT = 10, dN/dy = 1000) = 0.2. The RHIC 
blue curve is therefore independent of the b parameter. 
For larger a the blue curve flattens, while for smaller a 
the curve rises with pT more rapidly. The a = 1/3 rise is 
within the large error band of the PHENIX data shown 
in Fig. 1. 

It is important to note that if there were no density 
dependence of energy loss, then just because the spectral 
index decreases from ~ 5.5 to 4.5, the Raa{LHC) would 
be less quenched by - 70% at pr = 10 GeV/c at LHC! 
However, the approximate doubling of the initial density 
(and hence dE/dL) at LHC relative to RHIC results in 
halving Raa at pt = 10 GeV/c relative to RHIC. In 
addition, the reduction of the fractional energy loss with 
increasing pT with a = 1/3 causes Raa to rise with pT 
at LHC while at RHIC the px dependent increase of the 
spectral index compensates this natural rise. We will 
confirm these generic features in the next section with 
detailed WHDG numerical calculations. 

With this simple analytic model we can also easily 
explore the sensitivity of the LHC extrapolation to the 
path length dependence parameter b. The b parameter 
influences mainly the scaling exponent of dN/dy in the 
Bjorken expanding case. We see that decreasing b to 
(simulating perturbative elastic energy loss) decreases 
the difference between RHIC and LHC at px = 10 GeV/c 
but the change in px slope is similar. For b ~ 2 (simu- 
lating an path dependence suggested by AdS falling 
strings) predicts a significantly larger density dependence 
from RHIC to LHC. The cross over momentum where 
Raa at RHIC and LHC are equal increases monotoni- 
cally with b. The absolute value of Raa and its pr slope 
therefore provide strong constraints on the a, b parame- 
ters and therefore on the initial density or dN/ dy depen- 
dence, p(2-a-i-6)/3 jNjone of the parametric models predicts 
similar nuclear modification factors at RHIC and LHC 
above px > 10 GeV/c. 
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FIG. 2. Schematic nuclear modification factor model Eq. (8) 
for a "Bjorken plasma brick" normalized to Raa{pt = 10) = 
0.2 for central RHIC (blue) assuming a = 1/3. The momen- 
tum dependence at RHIC is independent of b once the refer- 
ence normalization at pr = 10 GeV/c is fixed. The red curves 
for LHC are evaluated with initial rapidity density scaled up 
by a factor 2.0 relative to RHIC. The 6 = 0, 1, 2 dependent 
extrapolations to LHC (red dashed, solid, dotted) are labeled 
by the equivalent L*' path length dependence of the total en- 
ergy loss in the static plasma limit. 



III. JET TOMOGRAPHY: QUANTITATIVE 

The extraction of tomographic (density dependent) 
information from Raa, and suppression observables in 
general, is complicated by many competing physics 
processes — some of which are interesting in their own 
right. Of greatest import for our discussion is the ini- 
tial phase space distribution of high-py quarks and glu- 
ons in both coordinate and momentum space. At RHIC, 
extensive reference spectra from p -I- p collisions at the 
same ^/s as in Au + Au collisions allow for a precise 
calibration of Raa and constrain experimentally the ini- 
tial quark and gluon spectra. The control d -I- Au data 
were essential to calibrate the magnitude of the initial 
state nuclear distortions of the initial jet spectra due to 
gluon saturation/shadowing physics [11, 61, 67, 83, 84]. 
At midrapidity the initial high-p^ partonic spectrum in 
the {x,Q'^) - (0.1, 100) - (0.2,400) kinematic range was 
not found to deviate significantly from the reference p -I- 
p pdfs (see Fig. 1 of [84]). Direct photon measurements 
further confirmed that the binary collision scaling based 
on Glauber diffuse nuclear reaction geometry can be used 
to understand the nuclear number and impact parameter 
dependence of hard processes [85]. 
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Comparable control measurements do not yet exist at 
LHC, thereby severely limiting the strength of conclu- 
sions that can be drawn from the first Raa data [f]. The 
absence of reference p + p data leads alone to about a 
factor ~ 2, p^-dependent systematic uncertainty in the 
normalization of the reported Raa- This uncertainty is 
shown for the 0-5% centrality data by the yellow band 
in Fig. 1. Second, the absence of control p + Pb as yet 
makes it impossible to deconvolute initial state nuclear 
suppression of high-p^ particles due to small x gluon 
saturation, or Color Glass Condensate (CGC) [11, 62- 
68], effects at LHC in the > 100 GeV^ range. Be- 
cause the fractional momenta relevant for midrapidity 
jet production at LHC are 10 times smaller at a given pT 
than at RHIC strong initial state suppression of the nu- 
clear/gluon structure has been predicted in [G8] at LHC. 

In this work, we assume the absence of significant 
initial-state suppression in the px ^ 10 — 20 GeV/c 
kinematic region corresponding to (x, Q^) ^ (0.01 — 
0.02, 100 - 400 GeV^). This is consistent with the 
DGLAP evolution of global fits to nuclear pdfs (see 
Fig.l of [84]). The first ATLAS measurement [86] of Z 
boson candidates is also consistent with unshadowed bi- 
nary scaling of jets in the x ~ 0.05 kinematic range. 
Future direct photon measurements at LHC will provide 
additional control over initial state shadowing/CGC ef- 
fects. 

The new ALICE data show features that appear strik- 
ingly similar to expectations based on pQCD energy loss: 
in particular Raa rises significantly at LHC as a function 
of pt rather than the observed fiatness within errors at 
RHIC. In [46] the predicted stronger rise as a function 
of momentum at LHC can be understood from the qual- 
itative model in Section II and the also from the follow- 
ing even simpler schematic model: The fractional energy 
loss of a high-pT- parton decreases in pQCD with mo- 
mentum as e ^ log{pT)/pT- The final momentum, p^, is 
related to the initial momentum, p^, by p;^ = (1 — £)Pt- 
For particle production distributions approximated by a 
power law, dN/dpT ^ Pr^ , the nuclear modification fac- 
tor Raa ^ ((1 — f)""^)- The suppression at RHIC is 
flater than at LHC due to an accidental cancellation be- 
tween 1) the fraction of high-pT gluons to quarks (see the 
upper panel of Fig. 3), 2) the hardening of the produc- 
tion spectrum as a function of pT (see the lower panel of 
Fig. 3), and 3) the decrease in energy loss as a function 
of Pt (see Fig. 4). 

As shown in the lower panel of Fig. 3 at LHC, the 
production spectrum is much harder (smaller ^ constant 
spectral index) than at RHIC. For a constant n{pT) the 
decrease of fractional energy loss with pt would lead to 
an Raa that increases with px- One can see from Fig. 1 
that the WHDG prediction rises with pr similar to the 
ALICE data. 

As one can see from the top panel of Fig. 3, also in con- 
trast to RHIC, LO pQCD predicts that hard jets at LHC 
are predominately gluons to much higher px- Naively 
then, with larger density, larger medium size (for Pb vs. 
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FIG. 3. (Top) Fraction of jets that are gluons as a function 
of jet PT from LO pQCD [87] for ^/s = 0.2, 2.76, and 5.5 TeV 
shown in red. Black curves show fraction of pions fragment- 
ing from gluons at final pion pr- (Bottom) Partonic spectral 
index, n{pT) ~ —d\og{dN / dydpT) / dlogijpT) , of initial y = Q 
gluon (red) and light quark (blue) jets compared at RHIC and 
LHC mp + p. 



Au), and the greater preponderance of gluons with 9/4 
enhanced energy loss relative to quark jets, should lead 
to an Raa suppressed much below that seen at RHIC. 
However, the smaller and fiatter spectral indices of both 
quarks and gluons compensates in the other direction. 

The first numerical GLV prediction for Raa i'^ 2002 
for 5.5 ATeV Pb -I- Pb collisions LHC including only ra- 
diative energy loss was given in Fig. 3 of [46] and over- 
laps well with the yellow ALICE systematic error band in 
Fig. 1. The first predicted WHDG Raa for LHC (see Fig. 
83 of [71]) including clastic energy loss as well as radia- 
tive energy loss and also realistic geometric jet path fluc- 
tuations accidentally remained close to Vitev's original 
GLV prediction. Our currently updated WHDG study 
incorporates the observed 2.2 times increase in charged 
particle rapidity density by ALICE in Pb-hPb 2.76 ATeV 
collisions to extrapolate the most recent 1 — ct uncertainty 
band of initial sQGP densities compared to high statis- 
tics PHENIX data from RHIC: this is the blue band in 
Fig. 1, which shows a significant underprediction of the 
ALICE data. 

We discuss in more detail below a range of theoreti- 
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cal and experimental issues that need further attention 
as jet tomography advances toward a more quantitative 
leveL The RHIC constrained WHDG results for Raa ap- 
pear consistent with the ALICE data and the earlier esti- 
mates within the present very large error band due to the 
unknown p + p baseline at this energy. While stronger 
conclusions are not warranted at this time, perturbative 
QCD tomography is at least consistent with the striking 
positive slope pr dependence of the nuclear modification 
light parton jet fragments at LHC. The change of the 
Pt slope, cIRaa/ dpT, from negative at SPS to approxi- 
mately zero at RHIC and finally to positive at LHC is a 
critical consistency test for the pQCD framework that jet 
holographic approaches must (eventually) also be able to 
pass. 



IV. FROM ENERGY LOSS TO Raa 
A. 0-5% Central 

The WHDG energy loss calculation described in [53, 
88] uses the first order in opacity radiative energy loss 
and Braaten-Thoma pQCD collisional partonic energy 
loss. The model assumes Debye-screened color scattering 
centers, and one loop screening mass /i = gT . Both light 
quarks and gluons have one loop medium-induced ther- 
mal masses of order /i. As noted previously, a generic 
feature of this pQCD radiative energy loss is that the 
fraction of final momentum to initial momentum of the 
parent parton, e = P^/p^' scales as e ^ \og{pT)/pT as 
can be seen Fig. 4. Of crucial importance in this formal- 
ism is the quantum interference between the hard vacuum 
radiation due to the initial creation of a hard parton and 
the radiation induced by scatterings in the medium. 

The redistribution of energy via collisional processes is 
found through the use of the fluctuation-dissipation theo- 
rem with the mean loss given by the thermal field theory 
calculation of Braaten and Thoma[89]. This redistribu- 
tion of energy can lead to the high-p^ parton cither los- 
ing or gaining energy, which is taken into account. The 
simple Gaussian approximation due to the fluctuation- 
dissipation theorem is a simplification that future more 
detailed but unwieldy elastic energy loss calculations can 
improve [90]. One can see in Fig. 4 that, with these as- 
sumptions and at medium densities appropriate for RHIC 
and LHC, at moderate px ^ 20 GeV elastic energy loss 
is as important as inelastic at both RHIC and LHC. 

We note in passing that the production points of the 
high momentum particles are distributed according to 
the binary distribution Taa] the medium density is as- 
sumed proportional to the participant density given by 
the Glauber model, and 1-D Bjorken-expansion is ap- 
proximately included by evaluating the density at a time 
one half the path length. 

Raa is defined as the ratio of observed spectra in A + 
A collisions for a given centrality divided by the observed 
spectrum in p p collisions scaled by the number of 
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FIG. 4. Average fractional energy loss (e) = (AE/E) for a 
light quark of initial momentum pr in a Bjorken-expanding 
brick of plasma of length 2 fm (blue) and 5 fm (red) for RHIC 
(top) and LHC (bottom) conditions from the WHDG energy 
loss model [53]. (e) is shown for purely radiative energy loss 
(Rad), purely collisional energy loss (Coll), and the incoherent 
convolution of the two (Conv). Elastic energy loss is a crucial 
contribution to total energy loss at both RHIC and LHC. 



binary collisions at the given centrality. ALICE measured 
the Raa of hadrons using the non-invariant spectra. 



I dprdy 



Ncoiiib)dNl^p+''' /dprdy 



(9) 



When one assumes that the production points of hard 
particles is proportional to the Taa distribution, where 
-^coU = / d^xT^p^lx), Ncow drops out of the theoretical 
calculation. However it is important to note that the 
values of iVcoii that we find using the optical Glauber 
model with the same parameters listed in [1] significantly 
disagree with those used by ALICE, found using a Monte- 
Carlo approach; we quantify this discrepancy in Table I. 
We will come back to this difference when we discuss R, 
below. 

In order to make contact with the experimen- 
tally observed hadrons, the partonic energy loss de- 
scribed above must be convolved with a partonic pro- 
duction spectrum and a set of fragmentation func- 
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ALICE 


ALICE 


Opt. Gl. 




Cent. 


A^coll 


Rel. Unc. 


A^coll 


% Diff 


0-5% 


1690 ± 131 


-8% 


1710 


~ 1% 


70-80% 


15.7 ±0.7 


-5% 


12.6 


-25% 



TABLE I. Values of A^coii used by ALICE and those found us- 
ing an optical Glauber model with a Woods-Saxon geometry 
and inelastic cross section identical to those used by ALICE in 
their Monte-Carlo calculation. Cf. to the uncertainties shown 
in Fig. 6. 



tions (FFs). We compute partonic suppressions via 
-^AA = ((1 - </)))"'■"'' (P^)-i), where we take = 
d\og{dN'i'^ / dpT) / d\og{pT) to approximate the power 
law for the partonic production spectrum, and (• • • ) rep- 
resents an averaging over the geometry for a given cen- 
trality class. Then 



^aa(pt; h) 



z dpT ^ z ' z ' 

dz dN'i ,pT.^^^h 
z dpT z 



■\z)+q^g/ 
(z) +q^g (10) 



We use a LO pQCD code [87] with a K factor of 3, 
CTEQ5L parton distribution functions [91] evaluated at 
Q = VtI'^-, and KKP fragmentation functions [92]. It is 
important to emphasize that the factor of two system- 
atic uncertainty in the K factor drops out of the ratio 
above. The produced neutral pion spectra expected for 
p -I- p collisions at RHIC and LHC energies using this LO 
pQCD procedure are shown in the top panel of Fig. 5. 
The LO pQCD spectrum at ^/s = 0.2 TeV agrees well 
with the PHENIX data to within ±20%. On the other 
hand we find that our LO pQCD spectrum with K = 3 
evaluated at y/s = 2.76 TeV systematically ovcrpredicts 
the spectrum used by ALICE by a nearly pT-indepcndent 
factor of ~ 2. We will come back to the absolute cross 
section normalization differences below. 

One of the key ingredients in an energy loss calculation 
is the density of the medium through which a fast parton 
propagates. We assume the transverse coordinate depen- 
dence of the density of the quark-gluon plasma medium 
is proportional to the participant density, dNg/d^xdy cx 
Ppart- The proportionality constant that connects these 
two quantities is precisely the lowest order tomographic 
information we can deduce by comparing energy loss 
calculations with data. The PHENIX collaboration ex- 
tracted this constant, reported as dNg/dy, via a rigorous 
statistical analysis comparing the WHDG energy loss to 
the 0-5% central RHIC RIa{pt) data [34]; it was found 
that the best fit value and one standard deviation un- 
certainty are dNg/dy = WOOtgy" for a fixed as = 0.3. 
Once this constant is fixed and we make the assump- 
tion that the quark-gluon plasma medium density scales 
with the observed number of charged hadrons, the QGP 
medium density at LHC is predicted. ALICE reported 
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FIG. 5. (Top) Hadron production spectrum data in p + p 
collisions from PHENIX at = 0.2 TeV [93] and from CMS 
at ,/s = 2.36 TeV [94]; the Tsallis distribution that is a best 
fit to the low pt CMS data is shown in red. Also in the figure 
is the LO pQCD calculation at ^/s — 0.2 TeV (black) and 
— 2.76 TeV (dashed black). The interpolation spectrum 
used by ALICE to normalize Raa at -^s = 2.76 TeV is shown 
in blue. (Bottom) The PHENIX data normalized to = 0.2 
TeV LO pQCD (black), the CMS Tsallis fit at = 2.36 
TeV normafized to LO pQCD at = 2.76 TeV (red), and 
the ALICE interpolation normalized to LO pQCD at = 
2.76 TeV (blue). Despite a larger y/s, note that the ALICE 
interpolation is a factor > 2 smaller than the CMS Tsallis fit, 
emphasizing the obvious need for 2.76 pp reference data. 



that the 0-5% central value of {dNch/drj}/ {Npiii-t/2) in- 
creased by a factor of 2.2 from RHIC [2]; the scaled den- 
sity in 0-5% central collisions at ^/snn 



2.76 TeV is 



+500 
-800- 



then dNg/dy = 3000 

With the QGP medium density at LHC so fixed, the 
suppression from the WHDG model is a constrained pre- 
diction; i.e. there are no free parameters in our cal- 
culation. We compare the resulting suppression for 0- 
5% central Rfp^ipT) in Fig. 1 to the 0-5% central AL- 
ICE charged hadron suppression. The results from the 
one standard deviation uncertainty in dNg/dy are rep- 
resented by thin blue lines; a light blue shaded region 
connects the two, denoting the uncertainty in the the- 



oretical calculation due to the uncertainty of the ex- 
tracted medium parameter from RHIC data. The dashed 
blue line represents R\j^{pt) for the best fit value of 
dNg/dy = 3000. One can see from the figure that the 
perturbative calculation qualitatively describes the in- 
crease in Raa as a function of pT, as we expected. Also 
our pocket formula from above describes the low-p^ nor- 
malization of the WHDG Raa results rather well. The 
small, < 0.1 value of Raa at low-p^ demonstrates again 
that the WHDG energy loss model is not fragile. The 
realistic distribution of production points and medium 
density means there is no (appreciable) corona of jets 
that emerge unmodified; there is no lower bound to the 
theoretical value of Raa- 

On the other hand, the very low normalization of the 
WHDG Raa seriously underpredicts the central values 
of the ALICE data. As there is no measured p p base- 
line at = 2.76 TeV at LHC yet, ALICE reports a very 
large uncertainty in Raa due to the uncertainty in the 
interpolated p -I- p baseline used. The result is that the 

0- 5% central WHDG energy loss calculation and the AL- 
ICE data agree at the edge of their respective reported 

1- (T uncertainties. It is worth emphasizing again, that 
while the WHDG Raa prediction is independent of the 
normalization of the p-l-p reference, the ALICE Raa is 
sensitive to the unmeasured p-|-p reference that has an in- 
trinsic theoretical factor of ^ two systematic uncertainty 
at both LO and NLO level. 

In addition to the theoretical uncertainty due to the 
finite precision extraction of a medium density at RHIC, 
there are also systematic uncertainties in the theoreti- 
cal calculation due to simplifying approximations made 
in the derivation of energy loss formulae. Some of these 
uncertainties have been quantified [95], and it turns out 
that the uncertainty due to the collinear approximation 
is in fact quite large: an exploration of this system- 
atic theoretical uncertainty yields a factor of 3 uncer- 
tainty in the extracted medium density at RHIC when 
energy loss is assumed to be radiative only. Given this 
large uncertainty, one might naively expect a similarly 
large uncertainty in the WHDG predictions for Raa at 
LHC. However this expectation is incorrect for two rea- 
sons: first, the uncertainty due to the collinear approx- 
imation decreases significantly when elastic energy loss 
is included [95]. Second, obscrvables such as V2{pt) — 
once the reference opacity is fixed — are found numeri- 
cally not to depend much on the specifics of the explored 
collinear approximation uncertainty [96] . The theoretical 
systematic uncertainties of the elastic energy loss contri- 
butions arc less well known and call for more scrutiny. 
In [90] it was shown that, contrary to often assumed 
multi soft Gaussian transverse clastic diffusion approx- 
imation, the large momentum transfer power law tails 
with q ~ lOr ~ 2 GeV momentum exchange (as included 
in WHDG through the \og{E/T) enhancement of elastic 
energy loss), cannot be ignored as assumed in many other 
models. In order to quantify the size of this collinear ra- 
diation approximation systematic uncertainty we use the 
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FIG. 6. Rep, the ratio of 0-5% central Raa(pt) to 70-80% pe- 
ripheral Raa{pt), from ALICE and our WHDG calculations. 
Experimental errors were added in quadrature, with those 
uncertainties due to the unknown p + p baseline cancelled. 



inelastic energy loss formula derived in Minkowski coordi- 
nates and take the maximum angle of emission for radia- 
tion to be 6'inax = 7r/2 (for a more detailed discussion, see 
[95]). The best fit value for the average medium density 
constrained by central RHIC data is dNg/dy = 1400 [95], 
which extrapolates to dNg/dy = 2800. The purple curves 
in Fig. 1 and denoted by "a; = xe" in Fig. 6 shows the 
values of RAAiPr) that result. That the curve lies within 
the WHDG uncertainty blue band due to the 1-cr medium 
parameter extraction at RHIC confirms that the system- 
atic theoretical uncertainty due to the collinear approxi- 
mation is small for Raa at LHC once the medium param- 
eter is fixed to RHIC data; in particular, the WHDG pre- 
dictions and the ALICE data are in quantitative tension 
within the combined experimental and currently quanti- 
tatively explored theoretical uncertainties. 



B. 70-80% Peripheral and R^p 

Continuing with the assumption that the QGP 
medium density scales with the number of charged 
hadrons we may make another constrained prediction 
from WHDG for the Raa of the 70-80% ccntrality 
class. From the ALICE centrality-dependent multiplic- 
ity data [4] we find that dNg/dy{70 - 80%) = 66^1°, 
and dNg/dy = 62 for the x = xe calculation. These are 
admittedly very small densities distributed over a small 
region, and it is not clear that a QGP medium forms in 
these highly peripheral heavy ion collisions. Neverthe- 
less we have assumed that our energy loss formalism is 
valid in both deconfined and confined media; it is not 
unreasonable to compare our results to the data, and our 
calculation and the ALICE results shown in Fig. 1. 

Since there is so much uncertainty in the p 4- p baseline 
spectrum, we find it useful to examine Rep, which is the 
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ratio of central Raa to peripheral Raa, in this case 



RcpiPr) 



RAAiPT;0-5%) 
Raa{pt\ 70-80%] 



(11) 



A 7t° WHDG LHC Extrapolation 
• ALICE 0-5% 
O ALICE 70-80% 



We plot both the experimental values and our calcula- 
tion in Fig. 6. Care was taken to propagate the relative 
systematic and statistical errors in quadrature. Some 
component of the systematic uncertainty shown in the 
original ALICE figures [1] is due to a pr-dependent un- 
certainty from the unknown p -|- p reference spectrum, 
which cancels in the Rep ratio. The ALICE paper [1] 
quotes the systematic uncertainties not due to the p + 
p spectrum as 5-7% and 8-10% in the central and pe- 
ripheral bins, respectively. As a conservative estimate 
we take the upper values and add (in quadrature) the 
relative uncertainty due to the A'con normalization; this 
procedure yields a systematic uncertainty of 15%, rep- 
resented by the gray bands in Fig. 6. Again the light 
blue band represents the theoretical uncertainty due to 
the extraction of the medium density at RHIC, with the 
dashed blue curve representing the best fit value, and 
the thin purple curve represents the result when using 
the X = xe calculation discussed above. 

The theoretical calculation reproduces the observed 
PT-dependence of Rep quite well. The theory results 
are systematically about 2 standard deviations of system- 
atic uncertainty below the data; however one should note 
that much of this systematic error is correlated and pt- 
indcpendcnt. and a smaller peripheral iVcoii, as suggested 
by the results from an optical Glauber calculation, would 
yield an experimental Rep suppression significantly closer 
to our prediction. Precise observations of direct photons 
and Z bosons should help reduce this possible extra un- 
certainty on the number of binary collisions in highly 
peripheral events. Additionally, that the optical Glauber 
results deviate so significantly from the Monte Carlo re- 
sults in the 70-80% centrality bin provides a quantitative 
feel for the importance of geometry fluctuations for these 
highly peripheral collisions; these possibly large geome- 
try fluctuations are not taken into account in the energy 
loss calculations presented here. Nevertheless, we find 
that this large discrepancy is a challenge to the pQCD 
paradigm assumed by the WHDG energy loss calculation. 

In Fig. 7 we present our predictions for R\a measured 
in Pb + Pb coUisions at 2.76 at LHC out to pr = 40 
GeV/c. Fig. 8 shows our predictions for central and 
peripheral suppression of neutral pions at ,/snn = 5.5 
TeV. Assuming that charged particle multiplicity scales 
as (sTVAr)*^'^^, as ALICE reported by fitting current world 
data [2], and assuming medium density scales with the 
observed charged particle multiplicity, we have that at 
^ysNN = 5.5 TeV the extrapolation from the RHIC ex- 
traction is dNg/dy = 3700t^Qooj the light blue band in 
the figure represents the predicted suppression based on 
these medium densities, with the dashed blue curve rep- 
resenting the most likely value. The thin purple curve 
represents the result when using the x = xe prescription 
and an extrapolation of dNg/dy — 3500. We note that 
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FIG. 7. Predictions of K^_^ for Pb + Pb collisions at 2.76 
ATeV from the WHDG energy loss model out to pr = 40 
GeV/c compared to current ALICE measurements of R')^JI 
[!]■ 



these constrained predictions for 5.5 ATeV are qualita- 
tively similar but differ in detail from our earlier 2007 pre- 
dictions that assumed a smaller medium density range, 
dNg/dy = 1700-2900 [72]. These updated WHDG con- 
strained predictions are qualitatively similar also to the 
predictions from 2002 GLV [46] that assumed smaller 
opacities at 5.5 ATeV but neglected the competing effects 
due to elastic and radiative energy loss as well as path 
length fluctuations included in WHDG. Similar pertuba- 
tive overquenching was also noted in other energy loss 
model approximations in [70, 97, 98]. As emphasized in 
Section II overquenching is a generic robust prediction of 
density dependent dE/dL models. 



V. CONCLUSIONS AND DISCUSSION 

In this paper we compared 7r° Raa(j>t) and Rcpipr) 
from the WHDG energy loss model [53] at ^/SNN = 2.76 
TeV and the recent ALICE data on charged hadron sup- 
pression at LHC [1]. The WHDG model includes both 
radiative and collisional channels and jet path length fluc- 
tuations. We found that at the momentum range cur- 
rently probed at LHC, collisional energy loss is as im- 
portant as radiative energy loss and cannot be neglected. 
Previous work [99] showed that even at top LHC ener- 
gies and pT ^ 250 GeV, elastic energy loss makes up a 
sizable fraction, ~ 25%, of the total energy loss in the 
HTL pQCD picture of the sQGP. The results we present 
assume that the HTL QGP medium density scales with 
the global charged particle multiplicity. Our results are 
true predictions based on rigorous statistical constraints 
of the medium density using WHDG at RHIC to the 
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from the WHDG energy loss model. 



QGP conditions at 2.76 central Pb+Pb as constrained 
by global dNch/dy from ALICE: there are no free pa- 
rameters. While unconstrained models without elastic 
energy and path length fluctuations can fit any Raa by 
adjusting the opacity, consistency with extensive RHIC 
jet quenching data leads in WHDG to a prediction of 
overquenching relative to the first ALICE Raa data. 

We find that our results show a qualitative agreement 
with the momentum dependence of the central ALICE 
data and a quantitative agreement with the normaliza- 
tion and pT-dependence of the peripheral ALICE results, 
within the reported uncertainties. Under our assumption 
of the scaling of the medium density, though, we find a 
quantitative disagreement with the normalization of the 
suppression for the reported central values of the ALICE 
data for 0-5% central Pb -I- Pb collisions. Unfortunately, 
the very large systematic uncertainty in the ALICE data 
due to the unmeasured base p -I- p reference precludes 
strong conclusions at this time. Our calculated Raa is 
free from this normalization ambiguity, but the experi- 
mental systematic error is dominated by it, and our pre- 
dictions for central Raa arc not in disagreement with 
the data within this large systematic uncertainty. Only 
a p -I- p run at ^/s = 2.76 TeV at LHC will resolve this 
uncertainty. 

In order to reduce systematic uncertainties, we there- 
fore compared Rep within the WHDG framework. An 
Rep analysis from ALICE was not released yet, and there- 
fore we calculated this quantity from the available results 
including propagating the reported uncertainties. 

The WHDG results for Rep also show qualitative agree- 
ment with the pr-dependent shape of the data, although 
the WHDG calculations seem to somewhat under-predict 
the slope for the central and peripheral Raa{pt) individ- 
ually. On the other hand, since our central Raa calcu- 
lations are more suppressed than the Raa data while 
the peripheral results compare favorably, our Rep results 
also over-suppress the Rep data. This disagreement is 



at about the 2-a level compared to the 15% systematic 
error we estimate from the ALICE data. Much of the 
15% systematic uncertainty is due to the A^coii normal- 
ization of Raa and is therefore pT-independent and cor- 
related, reducing the tension between the experimental 
results and the theoretical predictions. It is also worth 
noting again that the A^coU used by ALICE differs sig- 
nificantly from that given by an optical Glauber calcu- 
lation. Our results would compare more favorably with 
the data if the optical Glauber iVcon calculation would 
be used. Precise experimental measurements of high-pj^ 
probes that interact only weakly with the QGP plasma, 
such as direct photons or Z bosons, is the only way of 
experimentally constraining the A^coU normalization used 
in Raa and Rep analyses. Z bosons will give a much 
cleaner picture as photons arc created not only directly 
in the initial hard scatterings but also via photon brcm- 
sstrahlung [100-102]. 

As hard photons and Z bosons come from quark-anti- 
quark annihilation, a measurement of R\a ^^^^ 
than 1 could be interpreted as due to either a reduction 
in the expected number of binary collisions or the reduc- 
tion in the availability of anti-quarks due to gluon satu- 
ration. In particular there are predictions of a reduction 
by ~ 50% in the initial hard spectra at the momenta cur- 
rently observed at LHC due to gluon saturation effects 
[68]; this reduction would naturally decrease as a func- 
tion of pt as X increases. As such, the pT-dependence 
of the ALICE data may be due, at least in part, to ini- 
tial state effects. A precise measurement of R\a '^'^ ^aa 
showing the pT-independencc of these ratios would be 
required to demonstrate the lack of strong initial state 
shadowing in the x > 0.01 and > 100 GeV^ range 
relevant to the present ALICE data in the pT = 10 — 20 
GeV range (as expected from, e.g., [84]). p -I- A coUisions 
would provide a superb independent test of the possible 
influence of gluon saturation, similar to the crucial d -I- 
A collisions needed at RHIC to disentangle initial state 
from final state effects. 

Since there arc no adjustable parameters for us, the sig- 
nificant tension between our results and the ALICE data 
is a failure to simultaneously describe the normalizations 
of both the RHIC and LHC Raa{pt)- One possibility 
is the sQGP produced at LHC is in fact more transpar- 
ent than predicted by perturbative QCD tomographic 
models with medium densities that scale with observed 
particle rapidity densities. 

Theoretical possibilities that could contribute to the 
apparent transparency (decreased opacity) of the sQGP 
relative to the WHDG extrapolation from RHIC to LHC 
include 

1. Baryon anomaly [103-105] 

2. Gluon feedback [106] 

3. Gluon to quark jet conversion [101] 

4. Gluon self energy [48, 49] 
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5. Is the jet-medium coupling reduced at LHC: 
a,{LHC) < a,{RHIC) ? 

Item 1 can be resolved when identified tt, K,p, A, . . . , fl^ 
high-pT Raa becomes available. Item 2, in which the 
bremsstrahlung emitted gluons are kept track of and pro- 
duce observed hadrons. may be a partial explanation, 
but estimates so far have neglected the energy loss of 
the radiated gluons themselves. Radiated gluons with 
formation length less than the size of the medium could 
also be strongly quenched. A detailed centrality depen- 
dence of the di-hadron correlations may be able to re- 
solve such nonlinear gluon shower feedback mechanisms. 
Item 3 is possible channel that can effectively reduce the 
gluon jet color charge for asymmetric — 0.5| ~ 0.5) 
qq conversion. However, estimates[101] so far have ne- 
glected strong interference effects of medium and vac- 
uum induced amplitudes in finite size plasmas and high 
Pt octet color coherence of pair production. Jet flavor 
triggers that could discriminate light quark and gluon 
jets are needed to determine experimentally if this mech- 
anism is responsible for the apparent transparency of the 
sQGP at LHC. Item 4 involves HTL gluon dispersion 
effects on both the induced and hard initial radiation as- 
sociated with jet production. In [48, 49] it was shown 
that the Ter-Mikayelian and transition radiation effects 
reduce in general the effective dE/dL. Comparison of the 
centrality (or path length L) dependence of Raa may 
help untangle such dispersion effects. 

Item 5 refers to the possibility that the surprisingly 
transparent sQGP at RHIC could be due a reduction 
of the effective jet-medium coupling between RHIC and 
LHC. In the WHDG analysis here, we assumed a = 0.3 
is the same at RHIC and LHC and the Raa difference 
only reflects the increase of the initial sQGP density by 
2.2. In [70] an average path length approximation using 
GLV [42, 43] was used to show that approximately the 
same effective opacity L/A w 6 can fit both RHIC and 
LHC Raa- In the HTL approximation 1/A oc Tjc? up 
to slowly varying logarithmic corrections, so at L/ A con- 
stant implies a is reduced by a factor of 2.2^/^ — 1.3. 
However, in [70] ois =- 0.3 was assumed to be constant 
and L varied. In WHDG the effective L is completely 
fixed the distribution of path lengths. As l^Erad « 
and Ai5e/ cx the theoretical calculation of suppression 
depends very strongly on the specific value of taken. 
A fit, i.e. a postdiction to the LHC Raa data, can be 
achieved in WHDG by making the coupling a free pa- 
rameter and reducing it from the = 0.3 that we take 
in this analysis; we refrain from such uncontrolled fitting 
in this paper. A quantitative estimate of the effect of 
allowing the coupling to run, and therefore possibly be 
smaller at LHC than at RHIC, requires higher order the- 
oretical derivations that do not currently exist; in fact, 
even the qualitative result of such higher order effects 
are difficult to estimate as radiative energy loss calcula- 
tions always include soft exchanges between the leading 
parton and the medium particles that, in principle can 
only be handled by nonperturbative techniques. While 



one hopes that these higher order effects become small 
at higher leading parton momentum, there is always in 
heavy ion problems a temperature scale T which is the 
same order of magnitude as Aqcd- In particular, factor- 
ization has never been proven for energy loss calculations 
in heavy ion collisions. Could jet coupling to the sQGP 
at LHC be in fact more weakly coupled than at RHIC? 
From the near equality of bulk differentail elliptic flow, 
V2{pt < 2), the answer would appear to be no. How- 
ever, for short wavelengths pr > 10 GeV jet probes the 
effective jet medium coupling could in fact be smaller 
at LHC than at RHIC. The key observable to test this 
possibility may be the high pT elliptic and higher flow 
moments [75-78], V2{pt > 10 GeV/c). This observable 
remains a key stumbling block already at RHIC for all 
HTL/pQCD based models including WHDG that under- 
estimate ^2(5 GeV/c < Pt < 10 GeV/c) significantly. If 
V2{pt > 10 GeV/c) at LHC turns out to deviate less from 
WHDG — even when a{LHC) is reduced to account for 
the near identity of RHIC and LHC Raa — then a flrmer 
case could be made that the sQGP at LHC is indeed 
more transparent to jets than expected. 

In contrast to the above dynamical effects that could 
contribute to an apparent reduction in opacity at 
LHC relative to our WHDG expected growth L/X oc 
[dNchl dyY^^ , there are other dynamical effects neglected 
in our WHDG analysis that could contribute to an ap- 
parent enhancement of the sQGP opacity at LHC: 

1. High Color Glass Condensate [68] 

2. Dynamic magnetic scattering [80, 107, 108] 

3. AdS/CFT holography [22, 23, 26, 27, 81] 

Item 1 can be constrained via dedicated p-|-A and A-l-A 
direct gamma and Z boson control experiments. Cross 
correlating light and heavy quark flavor jet quenching 
systematics would provide quantitative insight into the 
potential influence of Items 2 and 3. 

There are several other possible sources of uncertainty 
that we did not address here. For instance one might ex- 
pect that the energy loss in confined matter would reflect 
the different properties of a hadronic medium as com- 
pared to a deconfined plasma of quarks and gluons [109]. 
Presumably, though, the cold matter energy loss would 
be smaller than hot, and this would lead to a greater dis- 
crepancy with the Rep data. Additionally, the ordering 
of length scales 1 //i <C \nfj] ^ L assumed in the DGLV 
energy loss derivations is violated for short length paths 
that may contribute more substantially to the hadrons 
that are ultimately observed at LHC as compared to 
RHIC due to the significantly more dense medium. Data 
from additional centrality classes will help clarify the pos- 
sible role of final state confined matter effects and length 
scale ordering dependencies. We also mentioned previ- 
ously that better calculations of the elastic energy loss of 
high-pT partons exist and can additionally be improved 
on. Aside from the open physics issues listed above. 
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there is a continuing need to improve numerical evalua- 
tion algorithms to remove simplifying numerical approx- 
imations used in WHDG and other tomographic models. 
Due to the above considerations, experimental measure- 
ments of observables out to very high pT^ for which we 
demand that theoretical calculations provide a consistent 
picture of both the mono- and di-jct data, will be crucial 
for furthering our understanding of the energy loss pro- 
cesses in experimentally accessible quark-gluon plasma, 
and hence crucially important for qualitatively and quan- 



titatively determining the properties of these plasmas. 
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